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Regular tesselations of space are characterized through their Schlafli symbols {p, ^, r}, where each 
cell has regular p-gonal sides, q meeting at each vertex, and r meeting on each edge. Regular 
tesselations with symbols {p, 3, 3} all satisfy Plateau's laws for equilibrium foams. For general p, 
however, these regular tesselations do not embed in Euclidean space, but require a uniform back- 
ground curvature. We study a class of regular foams on which, through conformal, stereographic 
projection to R*^ define irregular cells consistent with Plateau's laws. We analytically characterize a 
broad classes of bulk foam bubbles, and extend and explain recent observations on foam structure 
and shape distribution. Our approach also allows us to comment on foam stability by identifying a 
weak local maximum of A 2 /V at the maximally symmetric tetrahedral bubble that participates in 
T2 rearrangements. 



From beer to bread, metals to mousse, and shav- 
ing cream to souffles, foams are space-and mouth-filling 
structures which arise in a broad array of applied con- 
texts. Interest in foams has grown [1, 2] as they have been 
used to model many other systems: from copolymers and 
fuzzy colloids [s", to the Edwards' picture of granular 
matter [5] , and even to space-time and cosmology [gI, . 
Unfortunately, the details of a foam's bulk interior have 
long evaded both direct experimental probe [8] and ex- 
haustive theoretical treatment [9]. At issue is the topo- 
logical complexity of the network of foam bubble borders 
and the interaction of these borders with light, along with 
the disorder promoted by the incompatibility of Plateau's 
laws with space-filling regular polyhedra. Recall that 
Plateau posited El that three-dimensional dry foams 
would satisfy two rules: 1) each vertex terminates four 
edges, each an angle 6 = cos~^(— |) apart from the other, 
and 2) each edge is the border of three faces, and the di- 
hedral angle between these faces is 27r/3. In this letter 
we study the properties of a large class of bubbles that 
extend the notion of "isotropic Plateau polyhedra" [12 1 
or "average AT-hedra" [13] which we will call "average 
polyhedra" in the following. We determine both the area 
A and volume V of these cells and find that specific area, 
A^ /V is nearly constant for the class of dodecahedral 
and tetrahedral bubbles, with the exception of a weak 
local maximum around the symmetric tetrahedron. Our 
more general approach is a new step toward understand- 
ing the statistics of random foams 0], and the land- 
scape near the symmetric 4-sided bubble offers insight 
into foam evolution at finite temperature. The projec- 
tions that we make are all Plateau-law-abiding foams. 

The story in two dimensions is qualitatively different: 
there the decoration theorem p^] and the von Neumann- 
Mullins growth law [l^ provide a coherent and useful 
description of foams. In three-dimensions, geometrical 
frustration cuts our understanding at its core - indeed, 
only recently was the growth law generalized above two 
dimensions [17] and there are no regular tesselations of 
which satisfy Plateau's laws. Happily, it is well known 




FIG. 1: A simple flat (though non-conformal) projection of 
the 120-cell into R^. Darker lines are projected closer to the 
origin, lighter lines further. Note the dodecahedral nature of 
the cells and their foam-like coordination. 



that, in some cases, this type of frustration may be lifted 
by the introduction of a background curvature to the sys- 
tem [18, 19]. In order to progress along these lines in the 
case of foams, it is expedient to consider Plateau's laws 
in a language that allows for a simple connection to the 
curvature. That language is the geometric language of 
regular tesselations of space 0. To classify such tes- 
selations we use the Schlaffl symbols, {p, ^, r}, where p 
is the number of edges on each face of the cell, q is the 
number of these faces meeting at each vertex of one cell, 
and r is the number of cells or faces meeting around each 
edge. By way of example, the only such tesselation al- 
lowed in flat space is the cubic honeycomb. Since the 
cells are cubes, p = 4 and ^ = 3; four cubes fit together 
around an edge, and hence the Schlaffl symbol is {4, 3, 4}. 
The cubic honeycomb does not satisfy Plateau's rules as 
four faces meet at an edge. The requirement of dihedral 
angles of 27r/3 implies that r = 3. Since each vertex 
must terminate four edges, it follows that each cell must 
have three such edges and thus q = 3. Thus we see that 



any regular space tesselation of the form {p, 3, 3} neces- 
sarily represents a valid foam. If we divide each p-gonal 
face into p isoceles triangles, each will have one angle of 
27t/p and two angles |. On a uniformly curved space 
of constant sectional curvature this triangle, A, lies 
in a two-manifold of curvature K. The Gauss-Bonnet 
theorem relates the excess angle of the triangle to the 
integrated curvature: 

— + 2 X i cos-^ (--^ -7r= f KdA (1) 

In the flat limit K = and thus p ~ 5.1043; real foams 
are known to have '^5.1 sides per face on average (sl. [jlj 
as do the average polyhedra [l2|,ll3| . Of course, p must be 
a positive integer and the theory of polytopes [20'] limits 
us to p < 5 (we shun the frumious horocycle - {6, 3, 3}). 
It follows that the regular polytopes can be embedded 
in a three- space with K > 0, i.e. the three-sphere, S^. 
These "polychora" are the analogs of the two-dimensional 
tilings of which correspond to the platonic solids. We 
choose to study two tesselations of S^, the so-called 5- 
ceh, {3,3,3} and the 120-cen, {5,3,3} (Figured]). The 
120-cell is of interest because this regular poly tope re- 
quires the least amount of background curvature to real- 
ize among the family in question. Meanwhile, the 5-cell 
is worthwhile because it features tetrahedral bubbles like 
those seen at "T2" topological rearrangements, where a 
foam bubble shrinks until it disappears, leaving behind 
a four- fold vertex [9|. Related constructions have been 
used to show the uniqueness of /c-bubble clusters in d 
dimensions [22[. 

But how can we make contact with physical foams re- 
siding in actual, flat space? We observe that the con- 
straints of Plateau's laws that forced the form {p, 3, 3} 
are all constraints on angles, hence mappings from the 
3-sphere to that preserve angles will also preserve 
Plateau's laws. We choose here to work with the con- 
formal stereographic projection, long familiar to cartog- 
raphers [23] where X = x/{l — w)^Y = y/{l — w)^ and 
Z = z/{l — w) where X, F, Z are cartesian coordinates 
in the target R^, and x^y^z^w are the coordinates in 
C R^ with ^y^ ^ ^w'^ = 1. 

In two dimensions, this projection is the standard 
stereographic projection of to R^ and we may project 
the {p, 3} Platonic solids (the tetrahedron, cube, and 
dodecahedron) onto the plane. For concreteness, we 
will consider the {3, 3} tetrahedron shown in Fig. [21 
the conformal projection preserves the 27r/3 angle be- 
tween edges and so we will find a two-dimensional cell 
that satisfies Plateau's rules [2^. With one vertex 
at the pole, (0,0,1), the remaining set of vertices is 
q- — U [l _V2l -1\ (— Pi -1\ ((\ 2^2 

>^^V 3' 3 ' 3 '^ V 3' 3 ' S/'V'-'' 3 ' 3/J 

and they are all connected by arcs of great circles. We 
have the freedom to choose from which point we project 
S^, and each projection will give us two-dimensional cells 



2 




FIG. 2: (color online). The foam-like tetrahedral covering 
of the sphere and its Reuleaux Triangle image under stereo- 
graphic projection from a vertex. The full family of possible 
conformal distortions of this triangle is accessible by direct 
geometric calculation. 



that can fit in a two-dimensional foam. If, for instance, 
we project through (0,0, 1), then the top vertex is pro- 
jected to infinity and the remaining three vertices become 

^' = {(^' (0' 7f Since great cir- 

cles map to circles or lines, the edges of the tetrahedron 
become arcs of circles and, by symmetry, we have the fa- 
mous Reuleaux Triangle [25], familiar in two-dimensional 
foams as the ideal corner bubble during T2 topological 
rearrangements. Were we to project from another point, 
or equivalently, rotate the sphere and project through the 
pole, we would find asymmetric foams. Moreover, it is 
not just the central bubble that satisfies Plateau's laws; 
the entire projection is guaranteed to have all the correct 
angles for an equilibrium foam. The central cell on which 
we focus could exist in any equilibrium R^ foam. 

Working in R^, we can calculate both the perime- 
ter and area of a general triangular bubble in a two- 
dimensional foam. Here Plateau's rules are greatly sim- 
plified: edges must be arcs of circles and must meet at 
an angle of 27r/3. In order to parameterize an arbitrary 
triangle, we set two of the vertices at (—1,0) and (1,0), 
and allow the third to float as (x,?/), where x and y are 
both non- negative. It is tedious but straightforward to 
calculate both the perimeter and the area as a function 
of X and y. The exact formulae are complicated and the 
results are better characterized graphically in Figure [3l 

The machinery carries through in exactly the same 
way in the three dimensional case as well, allowing us 
to produce a continuous family of space foam bubbles. 
In principle, this technique can produce bubble families 
from each of the {3,3,3}, the {4,3,3}, and the {5,3,3} 
with 4-faced, 6-faced, and 12-faced bubbles, respectively. 
The class of transformations that belong to the confor- 
mal group of R^ is restricted to solid body rotations 
and translations, reflections, dilations, and spherical in- 
versions; the degrees of freedom available in the stere- 
ographic projection precisely account for each of these 
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FIG. 3: (color online). Contour plot of /V as a function 
of the triangle coordinate {x,y). The Reuleaux triangle is a 
local maximum, surrounded by three deep valleys. 



possibilities. Generically, parallel translation of the hy- 
persphere being projected relative to the target hyper- 
plane yields translations. Normal translation gives rise 
to dilations. Rotations and reflections of the hypersphere 
about hyperplanes that lie in the target hyperplane give 
the solid body transformations. The rest of the rotations 
give the spherical inversions. Since the shape of a bub- 
ble or bubble cluster is sensitive only to inversions from 
among these possible transformation types, we can ig- 
nore the translational degrees of freedom and focus only 
on the out-of-hyperplane rotations. 

Guided by experiment [8], we choose first to examine 
the most prevalent from among these bubble types, those 
with twelve faces. As in the two-dimensional example, we 
will have the freedom to choose the coordinate orienta- 
tion of the 120-cell inside and the ability to choose the 
base point of the stereographic projection; either one of 
these choices completely exhausts the degrees of freedom 
available, leaving the other to parameterize the resulting 
family of bubbles. We choose to pin the base point of the 
projection to the north pole of S^, (0, 0, 0, 1), and will ro- 
tate the 120-cell inside to generate different bubbles. 
In so doing it becomes transparent to identify the param- 
eterization of the foam bubble family with the solid body 
rotation taking place on the 120-cell pre-image. Said an- 
other way, we may unambiguously index our final bulk 
foam bubbles by elements of SO (4). Note that, due to 
the high symmetry of the pre-image 120-cell, it will suf- 
fice to consider the action of the projection on only one 
cell, since one cell may be carried into the position of any 
other by some element of SO (4). In analogy to the two- 
dimensional case, planar regions in S^ map to spherical 
caps in R^. 

If we view our parameter space as a higher dimen- 
sional rotation space, we will need to specify an unro- 
tated, reference state. We choose this state to coin- 



cide with the highest symmetry target bubble, the av- 
erage polyhedra Il3|, whose vertices coincide with 
the vertices of a regular dodecahedron and whose faces 
are all equivalent pentagonal spherical caps. This high- 
symmetry target bubble may be obtained by orienting 
the pre-image 120-cell so that its vertices lie at of the 

permutations of (±2, ±2, 0,0), (±^5, ±1, ±1, ±1), 
(±(^,±(^,±(^,±0-2), and (±(^^±(^-\±(^-\±(^-i) and 
of the even permutations of (±0^, ±1, 0), 

(±75, 0), and (±2, ±1, ±0-^) where = 

(l + y5)/2 is the Golden mean. 

With this characterization of the reference 120-cell in 
hand, we may begin to apply elements of SO (4) and 
project down to flat space. In analogy with the example 
calculation shown above for the projection of the tetra- 
hedron from S^ we obtain the equations of the spheres 
in defining the faces of the image dodecahedral bub- 
ble. Qualitatively, rotating the pre-image 120-cell around 
one of the w planes results in a bubble that has been 
stretched, with a vertex or edge pulled farther away from 
the others. This stretching also controls the sign of the 
face curvature for the bubbles; we consider relatively 
small deformations in order to preserve the positivity of 
the curvature for all twelve faces. 

Taking the appropriate face and cell intersections, we 
numerically determined the volume and surface area of 
the projected bubble. Perturbing away from the maxi- 
mally symmetric bubble by rotating around the w-z plane 
by an angle r, we calculate the non-dimensionalized com- 
bination A 2 / V. One would expect this quantity to dom- 
inate the bubble's free energy and, through the Boltz- 
mann factor, control the likelihood of finding a bubble 
of this shape in an arbitrary dry foam. Interestingly, 
however, we find that to within the errors of our calcula- 
tion /V is essentially constant with respect to r (see 
Figured]). This result extends the observations of Hilgen- 
feldt et al who found that for their idealized Isotropic 
Plateau Polyhedra the non-dimensional surface area is 
largely insensitive to the number of faces [12]. In ad- 
dition, the demonstration that even "stretching" foam 
bubbles does not change their reduced surface area-to- 
volume ratio serves to support and strengthen the recent 
results of Kraynik et ai., who found in simulations using 
Surface Evolver that the average reduced surface area- 
to-volume ratio of native bulk foam bubbles is constant 
[l3 |. The implication, then, is that the bubble shape dis- 
tribution in a bulk foam is dominated by the topology of 
the network of Plateau borders, with little contribution 
from bubble geometry. 

Pursuing the tetrahedral target bubbles generated by 
{3,3,3} is qualitatively similar to what was done with 
the 120-cell. The unrotated, reference state vertices 
for the 5-cell in S^ are : (0, 0, 0, 1), (^0, 0, , 

fo -l^fl fl^fl -l^fl -l^fl -l\ and 

\^^'Y6' 4Y3' 4y'\^2Y2' 2y 6^ 4.^3^ 4^' ^^^^ 
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FIG. 4: (color online). Dimensionless surface area to volume 
ratios for dodecahedral (lower curve) and tetrahedral bub- 
bles (upper curve) distorted by rotations in the w-z plane 
by T followed by sterographic projection. Two different ro- 
tational starting points are used for the tetrahedra, with 
no significant difference in /V. Note the near indepen- 
dence of both bubble types' ratios on the mapping parame- 
ter, T, with the exception of the maximum and valley near 
T = for the tetrahedron. Likewise, the distorted shapes 
are still highly constrained by Plateau's laws. For a sphere, 
A 2 /V — Ga/tF ~ 10.635 attains its minimum. 



(~^\/f ' ~\\[i'> ~\\[\'> ~i) • Pi'oceeding as before, we 
perturbed the maximally symmetric tetrahedral bubble 
by rotating an angle r about the w-z plane. As is the 
case for the dodecahedral bubbles, we find a largely flat 
dependence of A"^ /V on r (Figured]). Furthermore, in 
agreement with past work [14], we find that the values of 
the reduced surface area-to- volume ratio for the tetrahe- 
dral bubbles are essentially indistinguishable from those 
of the dodecahedral bubbles, differing by less than one 
percent on average. It is noteworthy, though, that the 
tetrahedra are uniformly higher in A^ /V ^ consistent with 
the expectation that the closer a cell is to achieving 13.7 
faces, the more efficiently it will enclose its volume. 

Our numerical results in three dimensions also suggest 
there is a local maximum of A 2 / V at the maximally sym- 
metric Reuleaux tetrahedron. Whether the existence of 
this local maximum is an impediment to T2 transitions at 
finite temperature is a natural question to pursue. Since 
this feature is preserved in two dimensions, where j A 
as shown in Figure [3l it may be possible to do both ex- 
periment and simulation to address this issue. 

We have generated a continuous class of dodecahe- 
dral and tetrahedral bubbles via conformal projection of 
known polychora and have shown that the specific area 
of these bubbles is roughly constant over our parameter 
space. Our methods have also allowed us to uncover a 
possible local maximum of A 2 /]/ at the T2-event bubbles 
in both two and three dimensions. It remains an inter- 
esting open question as to whether this approach can be 
generalized to study {p, 3, 3} polytopes for non-integer p 



along the lines of the average polyhedra [l2|, ll3| . 

It is a pleasure to acknowledge discussions with B.C. 
Chen, D.J. Durian, O.L. Halt, J. M. Kikkawa, R.B. Kus- 
ner, E. Matsumoto, and V. Vitelli. This work was sup- 
ported by NSF Grant DMR05-47320, a gift from L. J. 
Bernstein, and a gift from H. H. Coburn. 



[1] K. Feitosa, O. L. Halt, R. D. Kamien, and D. J. Durian, 
Europhys. Lett. 76, 683 (2006). 

B. Dollet and F. Graner, J. Fluid Mech. 585, 181 (2007). 
G.M. Grason, B.A. DiDonna, and R.D. Kamien, Phys. 
Rev. Lett. 91, 058304 (2003); G.M. Grason, Phys. Rep. 
433, 1 (2006). 

P. Ziherl and R.D. Kamien, Phys. Rev. Lett. 85, 3528 
(2000). 

S. F. Edwards and R. B. S. Oakeshott, Physica A, 
157(3), 1080 (1989). 

J.R .Gott, HI, M. Dickinson, and A.L. Melott, Ap. J. 
306, 341 (1986); C. Criado and N. Alamo, Class. Quan- 
tum Grav. 21 2241 (2004). 
J. A. Wheeler, Ann. Phys.2(6), 604 (1957). 

E. B. Matzke, American Journal of Botany, 33, 58 
(1945). 

D. Weaire and S. Hutzler, The Physics of Foams, Claren- 
don Press (Oxford, 1999). 

J.A.F. Plateau, Statique Experimentale et Theorique des 
Liquides Soumis aux Seules Forces Moleculaires (Gau- 
thier Villiard, Paris, 1873). 
J.E. Taylor, Ann. of Math. 103, 489 (1976). 
S. Hilgenfeldt, A.M. Kraynik, D.A. Reinelt, and J.M. 
Sullivan, Europhys. Lett. 67, 484 (2004). 
M.E. Glicksman, Phil. Mag. 85, 3 (2005). 
A. M. Kraynik, D. A. Reinelt, and F. van Swol, Phys. 
Rev. Lett. 93, 208301 (2004). 

F. Bolton and D. Weaire, Phil. Mag. B 63, 795 (1991). 
J. von Neumann, in Metal Interfaces, ed. by C. Her- 
ring, pp. 108110 (American Society for Metals, Cleve- 
land, 1952); W.W. Mulhns, in Metal Surfaces: Struc- 
ture, Energetics, and Kinetics ed. by W.D. Robertson 
and N.A. Gjostein, pp 17-66 (American Society for Met- 
als, Metals Park, Ohio, 1963). 

[17] R.D. MacPherson and D.J. Srolovitz, Nature 446, 1053 
(2007). 

[18] J. P. Sethna, D. C. Wright, and N. D. Mermin, Phys. 

Rev. Lett. 51, 467 (1983); J. P. Sethna, Phys. Rev. Lett. 

51, 2198 (1983). 
[19] C. D. Modes and R. D. Kamien, Phys. Rev. Lett. 99, 

235701 (2007); Phys. Rev. E 77, 041125 (2008). 
[20] H. S. M. Coxeter, Regular Polytopes, Methuen k Co. Ltd. 

(London, 1948). 
[21] R. Kusner, Proc. R. Soc. Lond. 439, 683 (1992). 
[22] A.M. Amilibia, Asian J. Math. 5, 23 (2001). 
[23] Ptolemy, Planisphaerium (ca. 150); see P. Kunitzsch and 

R. Lorch, Maslama's Notes on Ptolemy's Planisphaerium 

and Relatex Texts (Miinchen, Verlag der Bayer ischen 

Akademie der Wissenschaften, 1994). 
[24] A different but similar approach was studied in A. 

Mughal and D. Weaire, arXiv: 0801. 2474 vl (2008). 
[25] F. Reuleaux, The Kinematics of Machinery: Outlines of 

a Theory of Machines, Macmillan & Co. (London 1876). 



[2] 
[3] 



[4] 
[5] 

[6; 

[7] 
[8] 

[10; 



[11 
[12: 

[13; 

[14 

[is; 
[16; 



